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Abstract 

We provide predictions for charmonium and open charm production in nuclear col- 
lisions at SPS/FAIR energies within the framework of the statistical hadronization 
model. The increasing importance at lower energies of A c production is demon- 
strated and provides a challenge for future experiments. We also demonstrate that, 
because of the large charm quark mass and the different timescales for charm quark 
and charmed hadron production, possible modifications of charmed hadrons in the 
hot hadronic medium do not lead to measurable changes in cross sections for D- 
meson production. A possible influence of medium effects can be seen, however, 
in yields of charmonium. These effects are visible at all energies and results are 
presented for the energy range between charm threshold and RHIC energy. 



1 Introduction 



Charmonium production is considered, since the original proposal more than 20 years ago 
about its suppression in a Quark-Gluon Plasma (QGP) [lj, as an important probe to de- 
termine the degree of deconfinement reached in the fireball produced in ultra-relativistic 
nucleus-nucleus collisions. In a recent series of publications [2"f5P] we have demonstrated 
that, in the energy range from top SPS energy (y/s^N ~ 17 GeV) on, the data on J /if) 
and ij}' production in nucleus-nucleus collisions can be well described within the statis- 
tical hadronization model proposed in In particular the, at first glance surprizing, 
rapidity dependence of the nuclear modification factor R~aa observed by the PHENIX 
collaboration [6] is explained as due to the statistical hadronization of c and c pairs at 



the phase boundary between quark gluon plasma and hot hadronic matterQ We note 
that extrapolation to LHC energy of these results yields [3 4J a rather striking centrality 
dependence and, depending on the magnitude of the cc cross section in central Pb-Pb col- 
lisions, possibly even an enhancement (R^a > -0 °f ^/V production due to hadronization 
of thermalized, deconfmed and in general uncorrelated charm quarks. 

In the present publication we explore, for the first time, the lower energy range from 
near threshold (^/snn ~ 6 GeV) to top SPS energy. The lower part of this energy range 
can be investigated in the CBM experiment [7] at the future FAIR facility. One of the 
motivations for such studies was the expectation |7.8j to provide, by a measurement of 
D-meson production near threshold, information on their possible modification near the 
phase boundary . Here we demonstrate that, because of the relevant mass and time scales 
involved, medium effects on D-meson production are likely to be very small. Furthermore, 
because of the dominance of baryochemical potential (coupled with the charm neutrality 
condition) at low energies, it turns out to be important to measure in addition to D-meson 
production also the yield of charmed baryons to get a good measure of the total charm 
production cross section. 

In section 2 we will discuss the various time scales relevant for charm, charmonium, and 
open charm hadron production and discuss their relevance for the applicability of the 
statistical hadronization model as well as for the study of possible medium effects in the 
charm sector. Section 3 will provide a brief review of the statistical hadronization model. 
Our results on open charm and charmonium production from low beam energies on will 
be presented in section 4. In section 5 we will introduce various possible medium effects 
on open charm hadrons and study their influence on measurable quantities from FAIR to 
RHIC energies before concluding with a brief outlook in the last section. 



2 On relevant time scales and medium effects 



In the original scenario of J / ip suppression via Debye screening pQ it is assumed that the 
charmonia are rapidly formed in initial hard collisions but are subsequently destroyed in 
the QGP. While it is clear that the production of a (colored) charm quark pair takes 
place at time t c5 = l/(2m c ) < 0.1 fm, the formation time of charmonium involves color 
neutralization and the build-up of its wave function. The relevant time scale has been 
studied early on PfTU] and is of order 1 fm. Similar arguments also apply for the production 
time of charmed hadrons and we expect comparable time scales as for charmonium. 

We note that, at SPS energy where the 'melting scenario' was originally studied, this time 
is in the same range as the plasma formation time. At SPS and lower energies, charmonia 
can be formed in the pre-plasma phase and must be destroyed in the plasma if suppression 
by QGP is to take place. 



l J lib J I ib dNf}^ u /dy 

RjIa is defined as RjIX = — .^pp /d an d relates the charmonium yield in nucleus-nucleus 

coll J/tp' " 

collisions to that expected for a superposition of independent nucleon-nucleon collisions. In 
this expression, dNj/^/dy is the rapidity density of the J/tp yield integrated over transverse 
momentum and N co u is the number of binary collisions for a given centrality class. 



At the collider energies of RHIC and especially LHC the plasma formation time is likely to 
be much shorter (comparable to t c5 ). Furthermore, the number of charm quark pairs can 
exceed 10 per unit rapidity (central collisions at LHC). Initially, the 'collider' plasma will 
be hotter than Tp, the temperature above which screening takes place, and no charmonia 
will be formed at all in the QGP. It is our view that the charm quarks will be effec- 
tively thermalized leading to an uncorrelated pool of c and c quarks. Once the plasma 
temperature falls below To charmonia can be formed in principle, as well as destroyed, 
but as is indicated by the studies performed in [3], their formation rate is likely to be 
low. This finally leads to the notion, expressed explicitely in the statistical hadronization 
model, that all charmonia are produced by (re-) combination of charm quarks at the phase 
boundary. We would like to emphasize that, in this scenario, the particular value of To 
which is much discussed in the recent literature [TT][T2] . is not very important. Models 
that combine the 'melting scenario' with statistical hadronization have been proposed 
|13j . Alternatively, charmonium formation by coalescence in the plasma [mi5|H6|TTT] as 
well as within transport model approaches [TSyTQ"] has been considered. 

Another issue to be considered is the collision time t co a = 2R/'~f cm , where R is the radius 
of the (assumed equal mass) nuclei and 7 cm is the Lorentz 7 factor of each of the beams 
in the center-of-mass system. At SPS and lower energies, •y cm < 10 and t co u > 1 fm 
for a central Au-Au or Pb-Pb collision, so collision time, plasma formation time, and 
charmonium formation time are all of the same order [ID] . Furthermore, the maximum 
plasma temperature may not exceed Tp. In this situation the formed charmonia may be 
broken up by gluons and by the high energy nucleons still passing by from the collision. 
In this latter case, cold nuclear suppression needs to be carefully considered, as discussed, 
e.g., in [20]|2T| . However, we note in this context that the charm quarks resulting from such 
break-up processes eventually have to hadronize, which might again lead to charmonium 
production at the phase boundary if the charm quarks are kinetically thermalized, as is 
assumed in the statistical hadronization model [5"l3"] . Consequently, our calculations, in 
which both charmonium formation before QGP production and cold nuclear suppression 
are neglected, may somewhat underestimate the charmonium production yield at SPS 
energies and below. In that sense the below calculated medium effects are upper limits 
for energies close to threshold. 

At collider energies there will be yet another separation of time scales. At LHC energy, the 
momentum of a Pb nucleus is p cm =2.76 TeV per nucleon, leadingto 7 cm = 2940, hence 
tcoii < 5 ■ 10 -3 fm. Even "wee" partons with momentum fraction_| x w = 2.5 ■ 10~ 4 will 
pass by within a time t w = l/(xp cm ) < 0.3 fm, and will not destroy any charmonia since 
none exist at that time. We consequently expect that cold nuclear absorption will decrease 
from SPS to RHIC energy and should be negligible at LHC energy. First indications for 
this trend are visible in the PHENIX data [22]. 

Given the various time scales it becomes clear from the above discussion that the sta- 
tistical hadronization model should become a quantitative tool to describe charmonium 
and open charm production at collider energies without the explicit need to take account 
of any charmonium or open charm hadron formation before the QGP is developed and 
of cold nuclear absorption effects. We note in passing that the issue of shadowing or sat- 

2 We choose this value since it corresponds to a wee parton energy of approximately the binding 
energy of the J/ip meson. Smaller x values are hence not relevant for the present considerations. 



uration effects is of an entirely different nature: within the framework of the statistical 
hadronization model we need to know the rapidity density for open charm production 
in nucleus-nucleus collisions. Using this quantity, which of course contains shadowing or 
saturation effects, as input we can then provide cross sections for the production of all 
open and hidden charm hadrons. 

Finally we would like to discuss the effect of possible in-medium changes of charmed 
hadrons on their production cross section. We start the discussion by recalling that 

<y C c = ^{cr D + a Ac + (7g c + ...) + {a Vc + a J/4 , + a Xc + ...) (1) 



because of charm conservation. As shown in [3] , annihilation of charm quarks can be fully 
neglected. In the above equation, a D is the total cross section for the production of any 
D-meson. The cross section a c5 is governed by the mass of the charm quark m c m 1.3 
GeV |23j, which is much larger than any soft Quantum Chromodynamics (QCD) scale 
such as Aqcd- Therefore we expect no medium effects on this quantity. Such a separation 
of scales is not possible for strangeness production, and the situation there is not easily 
comparable. 

The much later formed D-mesons, or other charmed hadrons, may well change their mass 
in the hot medium. The results of various studies on in-medium modification of charmed 
hadrons masses [SJ24.25 26 27 28 29.30J are sometimes contradictory. Within a QCD sum 
rule model, Hayashigaki [27] predicts for D mesons a 50 MeV mass decrease at normal 
nuclear density, while for J ftp meson the shift is much smaller (5 MeV). Friman et al. [29] 
have concluded that the widths are little affected in dropping mass scenarios, while Tolos 
et al. [H] assert that only the widths could be affected by the nuclear medium but not 
the (pole) mass of D mesons. It is not clear whether the mass changes are different for 
particles and antiparticles, as advocated in [25], or identical because the vector potential 
may not matter if one assumes production at an early stage [28]. The "indirect" effect 
on J ftp production was also investigated [26,27. 18,29. 16J. A large mass shift of J/ip has 
been recently advocated [31] . We note that excellent fits of the common (non-charmed) 
hadrons to predictions of the thermal model have been obtained using vacuum masses 
[32.33 31. An attempt to use modified masses for the RHIC energy [53] has not produced 
a conclusive preference for any mass or width modifications of hadrons in medium. On 
the other hand, some evidence for possible mass modifications was presented in the chiral 
model of [36] . 

Whatever the medium effects may be, they can, because of the charm conservation equa- 
tion above, only lead to a redistribution of charm quarks. In particular, if the masses 
of D-mesons are lowered by the same amount at the phase boundary, this effect would 
practically not be visible in the D-meson cross section. Although the charm conserva- 
tion equation above is strictly correct only for the total cross section we expect within 
the framework of the statistical hadronization model, also little influence due to medium 
effects on distributions in rapidity and transverse momentum. This is due to the fact 
that the crucial input into our model is dN^ uAu / dy and there is no substantial D-meson 
rescattering after formation at the phase boundary. Modification of D-meson masses at 
the phase boundary will, however, influence the production rates for charmonia: after low- 
ering of their masses the D-mesons will eat away the charm quarks of the charmonia but 



since the D- mesons are much more abundant, their own yield will hardly change because 
of total charm conservation. 

Much of the above argument about medium effects is essentially model-independent and 
applies equally well at all energies. Here we will consider various types of scenarios for 
medium modifications and study their effect within the statistical hadronization frame- 
work in the energy range from charm threshold to collider energies. 



3 Reminder of ingredients and assumptions of the statistical hadronization 



The statistical hadronization model (SHM) |5,3J assumes that the charm quarks are 
produced in primary hard collisions and that their total number stays constant until 
hadronization. Another important factor is thermal equilibration in the QGP, at least 
near the critical temperature, T c . While data at RHIC energy [37] suggest charm equi- 
libration, at lower energies, where the initial densities and temperatures are lower, the 
assumption of equilibration can be questionable. In this exploratory study we have nev- 
ertheless assumed full thermalization. 

We focus on the energy range v / J/v r /v= 7-200 GeV and perform calculations for central 
collisions of heavy equal mass nuclei, corresponding to N part =350. We neglect charmonium 
production in the nuclear corona [3], since we focus in the following on central collisions, 
where such effects are small. 

In the following we briefly outline the calculation steps in our model [5|3] . The model 
has the following input parameters: i) charm production cross section in pp collisions; ii) 
characteristics at chemical freeze-out: temperature, T, baryochemical potential, and 
volume corresponding to one unit of rapidity Va v =i (our calculations are for midrapidity) . 
Since, in the end, our main results will be ratios of hadrons with charm quarks nomalized 
to the cc yield, the detailed magnitude of the open charm cross section and whether to 
use integrated yield or midrapidity yields is not crucial. 

The charm balance equation (5], which has to include canonical suppression factors [38] 
whenever the number of charm pairs is not much larger than 1, is used to determine a 
fugacity factor g c via: 



Here iV* r is the number of initially produced cc pairs and I n are modified Bessel func- 
tions. In the fireball of volume V the total number of open (iV*^ = n^V) and hidden 
(N^ = n^V) charm hadrons is computed from their grand-canonical densities n^. and 
n*g, respectively. This charm balance equation is the implementation within our model 
of the charm conservation constraint expressed in eq{H The densities of different parti- 
cle species in the grand canonical ensemble are calculated following the statistical model 
[32|33f34j . The balance equation (J2]) defines the fugacity parameter g c that accounts for 



model 





deviations of heavy quark multiplicity from the value that is expected in complete chem- 
ical equilibrium. The yield of charmonia of type j is obtained as: Nj = g^Nj h , while the 
yield of open charm hadrons is: N { = g c N- h Ii(g c Nj^)/ I (g c N^)- 



Fig. 1. Energy dependence of 
the charm production cross sec- 
tion in pp collisions. The NLO 
pQCD values [39] are compared 
to calculations using PYTHIA 
and to data in pA collisions, 
taken from ref. [40]. Our ex- 
trapolations for low energies are 
shown with continuous lines, for 
total and midrapidity (da" cg /dy) 
cross section. The open square 
is a midrapidity measurement in 
pp collisions [31] • The dashed 
line with dots indicates a param- 
eterization of the measured en- 
ergy dependence of the J/ ij) pro- 
duction cross section [42J . 



As no information on the charm production cross section is available for energies below 
yfs=lb GeV, we have to rely on extrapolation. The basis for this extrapolation is the 
energy dependence of the total charm production cross section calculated in ref. [39] for 
the CTEQ5M parton distribution functions in next-to- leading order (NLO), as shown in 
Fig. [U We have scaled these calculations to match the more recent values calculated at 
y / i=200 GeV in ref. [33]. We employ a threshold-based extrapolation using the following 
expression: 

Occ = k(l - ^r/V~s) a {V^r/^Ts) b (3) 

with fc=1.85 fib, yjs t hr=4:.5 GeV (calculated assuming a charm quark mass m c =1.3 GeV), 
a=4.3, and 6=-1.44. The parameters a, b, k were tuned to reproduce the low-energy part 
of the (scaled) NLO curve. The extrapolated curves for charm production cross section 
are shown with continuous lines in Fig. [TJ Also shown for comparison are calculations 
with PYTHIA [30] • To obtain the values at midrapidity we have extrapolated to lower 
energies the rapidity widths (FWHM) of the charm cross section known to be about 4 
units at RHIC [33] and about 2 units at SPS [33]- With these cross section values, the 
rapidity density of initially produced charm quark pairs strongly rises from 1.1- 10~ 3 to 1.7 
for the energy range 1 /s/v r /v= 7-200 GeV. We note that the so-obtained charm production 
cross section has an energy dependence similar to that measured for J/tp production, 
recently compiled and parametrized by the HERA-B collaboration [32]. For comparison, 
this is also shown in Fig. [IJ The extrapolation procedure for the low-energy part of the 
cross section obviously implies significant uncertainties. We emphasize, however, that the 
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most robust predictions of our model, i.e. the yields of charmed hadrons and charmonia 
relative to the initially produced cc pair yield are not influenced by the details of this 
extrapolation. 

For the studied energy range, T rises from 151 to 161 MeV from v / J/v r /v=7 to 12 GeV and 
stays constant for higher energies, while //& decreases from 434 to 22 MeV [34J . The volume 
Vaj/=i at midrapidity continuously rises 1341 from 760 to 2400 fm 3 . Due to the very small 
number of charm quarks at these low energies, the canonical suppression factor (Ii/Iq) 
is very large, but strongly decreases from about 1/30 to 1/1.2 for a /satjv=7 to 200 GeV. 
Correspondingly, the charm fugacity g c increases from 0.96 to 8.9. 



4 Energy dependence of charmed hadrons yield 

Our main results are presented in Fig. [2j The upper panel shows our predictions for 
the energy dependence of midrapidity yields for various charmed hadrons. Beyond the 
generally decreasing trend towards low energies for all yields one notices first a striking 
behavior of the production of A c baryons: their yield exhibits a weaker energy dependence 
than observed for other charmed hadrons. In our approach this is caused by the increase 
in baryochemical potential towards lower energies. A similar behavior is seen for the 
baryon. In detail, the production yields of D- mesons depend also on their quark content. 

The differing energy dependences of the yields of charmed hadrons are even more evident 
in the lower panel of Fig. [21 where we show the predicted yields normalized to the number 
of initially produced cc pairs. Except very near threshold, the J/ip production yield per 
cc pair exhibits a slow increase with increasing energy. This increase is a consequence of 
the quadratic term in the J/ifj yield equation discussed above. At LHC energy, the yield 
ratio J/ip/cc approaches 1% [3], scaling linearly with <r c5 ; for details see [45]. The if)' yield 
shows a similar energy dependence as the J/ip, except for our lowest energies, where the 
difference is due to the decrease of temperature (see above). We emphasize again that this 
model prediction, namely yields relative to cc pairs, is a robust result, as it is in the first 
order independent on the charm production cross section. Due to the expected similar 
temperature, the relative abundance of open charm hadrons at LHC is predicted [45] to 
be similar to that at RHIC energies. 



5 Effects of in-medium modification of charmed hadrons masses 

We consider two scenarios^ for a possible mass change Am of open charm hadrons con- 
taining light, u or d, quarks: i) a common decrease of 50 MeV for all charmed mesons 
and their antiparticles and a decrease of 100 MeV for the A c and S c baryons (50 MeV 



3 The scenarios are constructed by modification of the constituent quark masses of light (u and 
d) quarks in the charmed hadrons by fixed amounts. Reducing, for example, the light quark 
masses by 50 MeV will lower D-meson masses by 50 MeV and the A C (H C ) mass by 100 (50) 
MeV. 
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Fig. 2. Energy dependence of charmed hadron production at midrapidity. Upper panel: absolute 
yields, lower panel: yields relative to the number of cc pairs. Note the scale factors of 10 and 
100 for J/jJj and ip' mesons, respectively. 

decrease for H c ); ii) a decrease of 100 MeV for all charmed mesons and a 50 MeV increase 
for their antiparticles, with the same (scaled with the number of light quarks) scenario as 
in i) for the baryons. Scenario i) is more suited for an isospin-symmetric fireball produced 
in high-energy collisions and was used in [28J, while scenario ii) may be realized at low 
energies. In both scenarios, the masses of the D s mesons and of the charmonia are the 
vacuum masses. We also note that if one leaves all D-meson masses unchanged but allows 
their widths to increase, the resulting yields will increase by 11% (2.7%) for a width of 
100 MeV (50 MeV). If the in-medium widths exhibit tails towards low masses, as has 
been suggested by [8], to first order the effect on thermal densities is comparable with 
that from a decrease in the pole mass. 

The results for the two cases are presented in Fig. [3] as yields relative to the number 
of initially-produced cc pairs. As a result of the redistribution of the charm quarks over 
the various species, the relative yields of charmed hadrons may change. For example, in 
scenario i) the ratios of D-mesons are all close to those computed for vacuum masses 
(Fig. [2]), but the A c /D ratio is increased. Obvious are for scenario ii) the changes in the 
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Fig. 3. Energy dependence of the yield of charmed hadrons relative to the charm quark pair 
yield for two scenarios of the mass change (left panel for scenario i), right panel for scenario 
ii), see text). For the D mesons, the full and open symbols are for particles and antiparticles, 
respectively. Note the factors 10 and 100 for the J ftp and ip' mesons, respectively. 

relative abundances of the D and D mesons, as well as for the charmed baryons and also 
the relative production yields of and D~ mesons are very different. This difference 
is the result of the asymmetry in the mass shifts for particles and antiparticles. The 
change in yield of Df mesons occur as a consequence of the charm neutrality condition. 
Overall, however, charm conservation leads to rather small changes in the total yields. In 
contra-distinction, the effect of mass changes of charmed hadrons is very significant for 
charmonia, in particular the yields are more affected at low energies. 

In Fig. H] we demonstrate that the total open charm yield (sum over all charmed hadrons) 
exhibits essentially no change if one considers mass shifts, while the effect is large on 
charmonia. This is to be expected from eq. EJ as D-meson and A c -baryon masses are 
reduced, e.g., the charm fugacity g c is changed accordingly to conserve charm. Since the 
D-meson and A c -baryon yield varies linearly with g c we expect little change, while yields 
for charmonia vary strongly, since their yields are proportional to g\. To demonstrate this 
we plot, in Fig. HI the relative change of the yields with in-medium masses compared to 
the vacuum case. For this comparison, we have added a third case, namely considering 
that the mass change of charmed baryons is the same as for the mesons. 

Despite showing similar trends, in our model the reduction of the J / ip yield due to reduced 
in-medium masses of open charm hadrons has a different origin than that investigated in 
previous studies [26lf2TFl 8 29 ,16]. Dissociation of J/ip in a gas of 7r and p mesons [26] and 
effects originating from opening of decay channels of ip' and Xc states into DD [2"9lfT6] 
were considered in hadronic scenarios, while Zhang et al. [T8] have investigated the effect 
of mass changes both in the partonic (c quarks) and in the hadronic (D mesons) stage. In 
our model all hadrons with open and hidden charm are produced at chemical freeze-out. 
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Fig. 4. Energy dependence of the relative change in the production yield of open charm hadrons 
and of J ftp meson considering different scenarios for in-medium mass modifications (see text). 

Under the assumption that chemical freeze-out takes place at the phase boundary [46] , the 
medium effects could be due to the onset of chiral symmetry restoration or rescattering 
with the constituents of the medium. At present it cannot be ruled out that at lower 
energies the phase boundary does not coincide with chemical freeze-out. In this case, 
in-medium effects would be due to rescattering in the dense hadronic phase. 



6 Conclusions 

We have investigated charmonium production in the statistical hadronization model at 
lower energies. An interesting result is that the yield of charmed baryons (A c , H c ) relative 
to the total cc yield increases strongly with decreasing energy. Below ^/saw=10 GeV, the 
relative yield of A c exceeds that of any D meson except D Q , implying that an investiga- 
tion of open charm production at low energies needs to include careful measurements of 
charmed baryons, a difficult experimental task. The charmonium/open charm yield rises 
only slowly from energies near threshold to reach ~1% at LHC energy. Note that this ratio 
depends on the magnitude of the charm cross section, further underlining the importance 
to measure this quantity with precision. Our study is the first one addressing comprehen- 
sively the charm redistribution in principle and under various assumptions of in-medium 
masses of charmed hadrons. Because of a separation of time scales for charm quark and 
charmed hadron production, the overall D meson cross section is very little affected by 
in-medium mass changes, if charm conservation is taken into account. Measurable effects 
are predicted for the yields of charmonia. These effects are visible at all beam energies 
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• Am=-100MeV(D),_+5_0MeV(D) 
(A c ,£ c : -200 MeV, A ,£ c : +1 00 MeV) ; 

a Am=-50 MeV (A C ,Z C : -50 MeV) - 



and increase slightly towards threshold. 
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